Why INLA?

• Bayesian Hierarchical models
  • Many latent Gaussian variables
  • Few parameters
• Monte Carlo methods impractical
• Approximate posterior marginals using Laplace expansions
  • Clever algebra, end up with posterior in denominator
  • Taylor expand posterior about its mode
  • Results in a Gaussian approximation
  • Some numerical integration needed

Advantages and Disadvantages

• Advantages
  • Accurate approximation
  • Fast computation of many posterior marginals
• Disadvantages
  • Does not provide full joint posterior
  • Slow with >4-6 parameters

Derivations

From Blangiardo and Cameletti (2015), Section 4.7

Include derivations of \(\tilde{p}(\boldsymbol{\psi}|\mathbf{y})\), \(\tilde{p}(\theta_{i}|\boldsymbol{\psi},\mathbf{y})\), and \(\tilde{p}(\theta_{i}|\mathbf{y})\); and explanations of mode-finding and CCD.

Simple Normal Example

Blangiardo and Cameletti (2015) section 4.9

Spatial Point Pattern: Trees in a Rainforest